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Abstract 



We have compared the decay rates of Kl and K$ to 7r + 7r~ and n tt final 
states using a subset of the data from the KTeV experiment (E832) at Fer- 
milab. We find that the direct- CP- violation parameter Re(e'/e) is equal to 
(28.0 ± 3.0 (stat) ± 2.8 (syst)) x 1(T 4 . This result definitively establishes the 
existence of CP violation in a decay process. 
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The neutral K meson system has been the subject of much study since it was recognized 
that the two strangeness states (K°, K°) mix to produce short- and long-lived kaons (K$, 
K L ). The unexpected discovery of K L — > mi decays in 1964 |I| revealed that CP (charge- 
parity) symmetry is violated by the weak interaction, and it was soon understood that the 
dominant effect is an asymmetry in the K°-K° mixing, parametrized by e. Ever since, there 
has been great interest in determining whether CP violation also occurs in the K — > tttt 
decay process itself, an effect referred to as "direct" CP violation || and parametrized by 
e'. This would contribute differently to the rates of K L — > -n + n~ versus K L — ► 7r°7r° decays 
(relative to the corresponding K s decays), and thus would be observable as a nonzero value 
of 



Re(e'/e) » i 



T(K L -> U + TT-)/T(K S -> 7T+7T- 

Y(K L - 7r°7r )/r(K 5 



The standard Cabbibo-Kobayashi-Maskawa (CKM) model || can accommodate CP vi- 
olation in a natural way with a complex phase in the quark mixing matrix. The earliest 
standard-model calculations of Re(e'/e) which gave values of order 10~ 3 ~ 10~ 2 , were 
done before the top quark mass was known and before the importance of certain diagrams 
was appreciated. Modern calculations depend sensitively on input parameters and on the 
method used to estimate the hadronic matrix elements. Most recent estimates have tended 
toward values near or below 10~ 3 , for example (4.6 ±3.0) x 10~ 4 [f| and (8.5 ±5.9) x 1CT 4 ||; 
however, one group has estimated a larger range of values, (17±io) x 10 ~ 4 10] ■ Alternatively, 
a "superweak" interaction || could produce the observed CP-violating mixing but would 
give Re(e'/e) = 0. Therefore, a nonzero value of Re(e'/e) rules out the possibility that a 
superweak interaction is the sole source of CP violation. 

The two most precise past measurements of Re(e'/e) are in only fair agreement: the 
Fermilab E731 experiment reported Re(e'/e) = (7.4 ± 5.9) x 1(T 4 f|, while CERN NA31 
found Re(e'/e) = (23±6.5) x 10~ 4 [|It]]. Because of the importance of definitively establishing 
the existence of direct CP violation and determining its magnitude, new experiments have 
been undertaken at Fermilab, CERN, and Frascati to measure Re(e'/e) with precisions of 
(1 ~ 2) x 1(T 4 . 

This Letter reports a new measurement of Re(e'/e) using 23% of the data collected by 
the KTeV experiment (E832) during the 1996-97 Fermilab fixed-target run. The KTeV 
experiment was designed to improve upon the previous generation of experiments and ulti- 
mately to have the sensitivity to establish direct CP violation in the range of the smaller 



estimates in || and ||. The experimental technique is the same as in E731 |TTJ, and differs 
from NA31 in two key ways. First, it uses two kaon beams from a single target to enable 
the simultaneous collection of Kl and Ks decays in order to be insensitive to the inevitable 
time variation of beam characteristics and detector inefficiencies. Second, it uses a precision 
magnetic spectrometer to minimize backgrounds in the 7r + 7r~ samples and to allow in situ 
calibration of the calorimeter. While the method of producing the Ks beam (by passing a 
Kl beam through a "regenerator") is also the same as in E731, the KTeV regenerator is 
made of scintillator and is fully instrumented to reduce the scattered-kaon background to 
the coherently-regenerated signal. A new beamline was constructed for KTeV with much 



cleaner beam collimation and improved muon sweeping [ 12j| . Finally, the KTeV electromag- 
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FIG. 1. Plan view of the KTeV apparatus as configured to measure Re(e'/e). The evacuated 
decay volume ends with a thin vacuum window at Z = 159 m. The label "Csl" indicates the 
electromagnetic calorimeter. 



netic calorimeter has much higher precision than the E731 calorimeter, permitting more 
accurate 7r°7r° reconstruction and better background suppression. 

Figure [1] shows the two beams (called "regenerator" and "vacuum") in the evacuated 
decay volume, with the main detector elements located downstream. The regenerator alter- 
nates sides between accelerator extractions to minimize the effect of any left-right beam or 
detector asymmetry. A "movable absorber," far upstream, attenuates the beam incident on 
the regenerator. To measure the double ratio of decay rates in the expression for Re(e'/e), 
we must understand the difference between the acceptances for K s versus K L decays to 
each 7T7T final state. Triggering, reconstruction and event selection are done with identical 
criteria for decays in either beam, so the principal difference between the Ks and Kl data 
samples is in the decay vertex distributions, shown in Fig. ||] as a function of Z, the distance 
from the kaon production target. Therefore, the most crucial requirement for measuring 
Re(e'/e) with this technique is a precise understanding of the Z-dependence of the detector 



acceptance (T^j 



In the regenerator beam, the beginning of the decay region is sharply defined by a 
lead-scintillator module at the downstream end of the regenerator. In the vacuum beam, 
the acceptance for decays upstream of Z = 122 m is limited by the "mask anti" (MA), a 
lead-scintillator counter with two square holes 50% larger than the beams. 

The KTeV spectrometer consists of four rectangular drift chambers, each with two hor- 
izontal and two vertical planes of sense wires, and a large dipole magnet which imparts a 
transverse momentum of 0.412 GeV/c. The spaces between the drift chambers are filled 
with helium to reduce scattering. The drift chambers measure horizontal and vertical track 
position with a typical resolution of 110 (ttrn and momentum with a resolution of 0.4% at 
the mean pion momentum of 36 GeV/c. 

The electromagnetic calorimeter [TJ|] consists of 3100 blocks of pure cesium iodide (Csl) 
in a square array 1.9 m on a side and 0.5 m deep. Two 15 cm square beam holes allow 
passage of the neutral beams through the calorimeter. The calorimeter was calibrated using 
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FIG. 2. Decay vertex distributions for the (a) K 




and (b) K 



decay modes, 



showing the difference between the "regenerator" (Ks) and "vacuum" (Kl) beams. 



1.9 x 10 8 momentum-analyzed electrons from Kl — > Trez/ decays collected during normal 
running. Individually tuned wrapping of the Csl blocks and the development of very-low- 
noise readout electronics |15| have enabled the calorimeter to achieve an average energy 
resolution of 0.7% for photons from 7r 7r° decays, which have a mean energy of 19 GeV. 

The inner aperture for photons at the Csl is sharply defined by a tungsten-scintillator 
"collar anti" (CA) counter around each beam hole. In addition, there are ten lead-scintillator 
"photon veto" counters to detect particles escaping from the decay volume or missing the 
Csl, in order to suppress Kl — > 3n background in the 7r 7r° samples. 

The trigger system initiates detector readout based on synchronous signals from a scin- 
tillator hodoscope located upstream of the calorimeter (for 7r + 7r~) or on a fast analog energy 
sum from the calorimeter (for 7r°7r°). To keep the trigger rate at a manageable level, triggers 
are inhibited by fast veto signals from the regenerator, the MA, a subset of the photon vetoes, 
and a downstream hodoscope located behind 4 m of steel to detect muons [[L6[]. For 7r + 7r~, 



additional requirements are made on the number and pattern of hits in the drift chambers. 
For 7r°7r°, a hardware processor |I7fl must find 4 or 5 "clusters" of energy in the calorimeter. 
After readout, a CPU-based "Level 3 filter" reconstructs events and applies some loose kine- 
matic cuts to select 7r + 7r~ and 7r 7r° candidates. Besides these signal modes, large samples 
of Kl — > 7rez/, Kl — > 7r + 7r~7r°, and Kl — > 3tt° decays are recorded for detector calibration 
and acceptance studies. In addition, an "accidental" trigger is formed, using scintillation 
counters near the kaon production target, to randomly record the underlying activity in the 
KTeV detector with the same instantaneous-intensity distribution as the physics data. 

The 7r°7r° samples used for this analysis are from the data collected in 1996, while the 
7r + 7r~ samples are from the first 18 days of data collected in 1997. We decided not to use 
the 7r + 7r~ data from 1996 because the Level 3 filter had a 22% inefficiency arising from 
an unanticipated drift chamber effect which sometimes delayed a hit by 20 ns or more. 
The inefficiency was nearly the same for both beams but would still have led to a large 
systematic error on Re(e'/e). The Level 3 software was modified for the 1997 run to allow 
for this effect, resulting in an inefficiency of less than 0.1%. Using 7t + tt~ and data from 
different running periods does not significantly increase the systematic error on Re(e'/e) 
because the two modes use essentially independent detector systems; detector inefficiencies 
and sources of deadtime cancel in the Ks/Kl ratio for either mode independently. The only 
direct effect is a possible difference in the Ks/K L flux ratio, which will be discussed later. 
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For offline selection of tt + tt~ candidates, each pion is required to have a momentum of 
at least 8 GeV/c and to deposit less than 85% of its energy in the calorimeter. In order 
to cleanly define the acceptance and to avoid topologies with poor reconstruction efficiency, 
cuts are made on the distance from each pion to the edges of the drift chambers, calorimeter, 
MA, and CA, and on the separation distance between the two pions at the drift chambers 
and calorimeter. The ir + 7T~ invariant mass is required to be between 488 and 508 MeV/c 2 
(where the mean resolution is approximately 1.6 MeV/c 2 ) and the square of the transverse 
momentum of the 7r + 7r~ system relative to the initial kaon trajectory, p 2 -, is required to be 
less than 250 MeV 2 /c 2 . 

After applying various corrections to the raw calorimeter information, 7r°7r° candidates 
are reconstructed from four-photon events by choosing the photon pairing combination which 
is most consistent with the hypothesis of two ir° decays at a common point, interpreted as 
the kaon decay vertex. Each photon is required to have an energy of at least 3 GeV and to 
be at least 5 cm from the outer edge of the Csl and 7.5 cm from any other photon. The 
four-photon invariant mass is required to be between 490 and 505 MeV/c 2 , where the mean 
mass resolution is approximately 1.5 MeV/c 2 . The initial kaon trajectory is unknown, so 
the only available indicator of kaon scattering is the position of the energy centroid of the 
four photons at the Csl. This is used to calculate a "ring number" , defined as four times the 
square of the larger normal distance (horizontal or vertical), in centimeters, from the energy 
centroid to the center of the closer beam. Its value is required to be less than 110, which 
selects events with energy centroid lying within a square region of area 110 cm 2 centered on 
each beam. 

In both the 7t + tt~ and tt°tt° analyses, cuts are made on energy deposits in the MA, photon 
veto counters, and regenerator. The final samples consist of events with 110 < Z < 158 m 
and 40 < E K < 160 GeV. 

A detailed Monte Carlo (MC) simulation is used to determine the detector acceptance for 
the 7T7T signal modes and to evaluate backgrounds. The simulation models kaon production 
and regeneration to generate decays with the same energy and Z distributions as the data. 
The decay products are traced through the KTeV detector, allowing for electromagnetic 
interactions with beamline material and for pion decay. The acceptance is largely determined 
by the geometry of the detector and by geometric analysis cuts; however, to understand 
reconstruction biases it is important to simulate the detector response accurately. Energy 
deposits in the Csl blocks from photons, pions, and electrons are based on the GEANT 
package [Tj|. Drift chamber inefficiencies and the delayed-hit effect are simulated using 



parametrizations and position dependences measured from tc + it~ data, ireu and 3tt° data 
samples are used to check or tune various aspects of the detector geometry and simulation. 
To reproduce possible biases due to underlying activity in the detector, an event from the 
accidental trigger is overlaid on top of each simulated decay; the net effect on the measured 
value of Re(e' / e) is of order 10~ 4 . MC event samples are subjected to the same reconstruction 
and selection criteria as the data samples. 

Background contributions to the 7r + 7r _ samples are determined by using sidebands in 
the mass and p\ distributions to normalize MC predictions from the various background 
processes. Figure [| (a) and (b) show that the "p\ distributions for data are well described 
by the sum of coherent nit MC and total background MC. Kl — > ireu and Kl — > t\\iv 
decays, with the electron or muon misidentified as a pion, contribute 0.069% (0.003%) to 
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FIG. 3. Distributions of pj> for the 7r + 7r _ samples and ring number for the 7r 7r° samples. 
Total background levels and uncertainties (dominated by systematics) are given for the samples 
passing the analysis cuts (arrows). 



the vacuum (regenerator) beam. The dominant regenerator-beam background (0.072%) is 
from kaons which scatter in the regenerator before decaying to 7r + 7r~. Kaons which scatter 
in the final beam-defining collimator contribute an additional 0.014% to each beam. Data 
samples of tc + tt~ decays from kaons which scatter in the regenerator or collimator are used 
to tune physics-motivated scattering models incorporated into the MC simulation. 

The background levels are much larger for the samples since the ring-number vari- 
able is not as effective as p\ at identifying scattered kaons and cannot detect "crossover" 
scattering from the regenerator into the vacuum beam. Ring-number distributions are shown 
in Fig. [| (c) and (d). The upturn under the peak in (c) is due to K L — > 3n° decays with 
lost and/or overlapping photons; it is determined, using mass sidebands, to contribute a 
background of 0.27% (0.01%) to the vacuum (regenerator) beam. A ring-number sideband 
(286-792) is used to normalize MC distributions from kaons that scatter before decaying to 
7r°7T . The vacuum (regenerator) beam background includes 0.30% (1.07%) from regenerator 
scattering and 0.16% (0.14%) from collimator scattering. Pairs of 7r°'s produced by hadronic 
interactions in the regenerator contribute an additional background of 0.01% in that beam. 

After background subtraction, the net yields are 2,607,274 tt + tc~ in the vacuum beam, 
4,515,928 7r + 7r _ in the regenerator beam, 862,254 7r°7r° in the vacuum beam, and 1,433,923 
7r°7r° in the regenerator beam. 

Re(e'/e) is extracted from the background-subtracted data using a fitting program which 
calculates decay vertex distributions, properly treating regeneration and Ks-K L interference 
(including the residual K$ component in the vacuum beam at high energy). The acceptance 
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TABLE I. Systematic uncertainties on Re(e'/e). 





Uncertainty (xlO 4 ) 


Source of uncertainty 


from 7r + 7r~ from 7r°7r° 


Class 1: Data collection 




Trigger and Level 3 filter 


0.5 0.3 



Class 2: Event reconstruction, selection, backgrounds 

Energy scale 

Calorimeter nonlinearity 

Detector calibration, alignment 

Analysis cut variations 

Background subtraction 
Class 3: Detector acceptance 

Limiting apertures 

Detector resolution 

Drift chamber simulation 

Z dependence of acceptance 

Monte Carlo statistics 
Class 4: Kaon flux and physics parameters 

Regenerator-beam attenuation: 



1996 versus 1997 0.2 

Energy dependence 0.2 

Am, ts, regeneration phase 0.2 

TOTAL 2^8 



0.1 

0.3 
0.6 
0.2 



0.7 
0.6 
0.4 
0.8 
0.8 



0.3 
0.4 
0.6 
1.6 
0.5 



0.5 
<0.1 



0.7 
0.9 



correction (as determined from MC) is applied, and the resulting prediction for each decay 
mode is integrated over Z and compared to the data in 10-GeV bins of kaon energy. CPT 
symmetry is assumed, and the values of the K^-Ks mass difference (Am) and K$ lifetime 
(ts) are fixed to the published average values [TE|. The regeneration amplitude is floated 
in the fit, but constrained to have a power-law dependence on kaon energy, with the phase 
determined by analyticity fTTJ^^l • The kaon energy distributions are allowed to be different 



for the tt + it~ and 7r°7r° modes, with a floating normalization correction in each energy bin for 
each mode (24 fit parameters). Fitting was done "blind", by hiding the value of Re(e'/e) with 
an unknown offset, until after the analysis and systematic error evaluation were finalized. 
The result is Re(e'/e) = (28.0 ±3.0) x 10~ 4 , where the error is statistical only. The fit x 2 is 
30 for 21 degrees of freedom. 

As a general rule, only biases which affect the Ki and K$ samples differently will lead 
to systematic errors on Re(e'/e). Possible sources may be divided into four classes: (1) data 
collection inefficiencies; (2) biases in event reconstruction, sample selection, and background 
subtraction; (3) misunderstanding of the detector acceptance; and (4) uncertainties in kaon 
flux and physics parameters. Table H summarizes all of the estimated contributions; only 
those that are large or require special explanation will be discussed below. 

Two of the largest uncertainties in the second class are related to the measurement 
of photon energies by the calorimeter. A systematic shift in measured energies can shift 
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the reconstructed Z vertex and Ek distributions for the 7r°7r° sample and thus can bias 
Re(e'/e), mainly by moving Ki events past the fiducial Z cut at 158 m. After calibrating 
the calorimeter with electrons (and allowing for a small expected electron-photon differ- 
ence), a final energy scale correction for photons of —0.125% is determined by matching the 
sharp turn-on of the distribution at the regenerator edge between data and MC. 

After making this correction, a check using tt° pairs produced by hadronic interactions in 
the vacuum window reveals a Z mismatch of 2 cm at the downstream end of the decay 
region, leading to a systematic error of 0.7 x 10~ 4 on Re(e'/e). Residual nonlinearities in 
the calorimeter response, studied from the variation of the mean 7r°7r° invariant mass as a 
function of Ek, contribute an additional error of 0.6 x 10~ 4 . 

We assign systematic errors based on the dependence of the measured value of Re(e'/e) 
on variations of key analysis cuts, in particular the p\ cut for ix + -k~ and the ring-number 
and photon quality cuts for 7r°7r°. No significant dependence is observed on other analysis 
cuts. 

The accuracy of the background determination for the 7r°7r° samples depends on our 
understanding of kaon scattering in the regenerator and collimator. We consider several 
variations in scattering models and in the procedures for tuning the MC with K — > 7t + tt~ 
decays; these affect the shapes of the background MC ring-number distributions, but the 
sideband normalization procedure limits the impact on Re(e'/e). We assign an uncertainty 
of 0.8 x 10~ 4 . 

The third class of systematic uncertainties, related to detector acceptance, contributes 
the most to the total systematic error. Many potential detector modeling problems would 
affect the acceptance as a function of Z, so a crucial check of our understanding of the 
acceptance is to compare the Z distribution for the data against the MC simulation. Figure |] 
shows the vacuum-beam comparisons for the 7t + tt~ and 7r°7r° signal modes as well as for the 
much larger Treu and 3tc° samples. The overall agreement is very good, but since the mean 
Z positions for Kl and K~s decays differ by about 6 m, a relative slope of 10~ 4 per meter 
in the data/MC ratio would cause an error of 10~ 4 on Re(e'/e). As shown in Fig. |] (b), the 
Treu comparison agrees to better than this level; however, the 7r + 7r~ comparison has a slope 
of ( — 1.60 ± 0.63) x 10~ 4 per meter. Although the statistical significance of the 7r + 7r~ slope 
is marginal, we assign a systematic error on Re(e'/e) based on the full size of the apparent 
slope, 1.6 x 10~ 4 . The 3n° and 7r°7r° Z distributions agree well, and we place a limit of 
0.7 x 10~ 4 on the possible Re(e'/e) bias from the neutral-mode acceptance. 

Other checks on the acceptance include data/MC comparisons of track illuminations at 
the drift chambers and Csl, photon illumination at the Csl, and minimum photon separation 
distance. These all agree well and indicate no other sources of acceptance misunderstanding. 

The final class of systematic uncertainties includes possible differences in the Ks/Kl 
flux ratio between the 7r + 7r~ and 7r 7r° samples. The flux ratio is nominally the same for the 
1996 and 1997 running periods because the same regenerator and movable absorber were 
used; however, we must assign a small uncertainty on Re(e'/e) due to a possible temperature 
difference which would change their densities and thus the regenerator-beam attenuation. 
In addition, the 7r + 7r~ and 7r 7r° samples have somewhat different energy distributions, so 
the uncertainty in the energy dependence of the attenuation (measured using tt + tt~7t° and 
37T° data) leads to a small uncertainty on Re(e'/e). 

Finally, we assign uncertainties corresponding to one-sigma variations of Am and t$ from 
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the published averages pi] , and from a deviation of the phase of the regeneration amplitude 



by ±0.5° from the value given by analyticity [20]. Adding all contributions in quadrature, 
the total systematic uncertainty on Re(e'/e) is 2.8 x 10~ 4 . 

We have performed several cross-checks on the Re(e'/e) result. Consistent values are 
obtained at all kaon energies, and there is no significant variation as a function of time or 
beam intensity. Relaxing the power-law constraint on the regeneration amplitude yields a 
consistent value with the same precision. We have also extracted Re(e'/e) using an alter- 
native fitting technique which compares the vacuum- and regenerator-beam Z distributions 
directly, eliminating the need for a Monte Carlo simulation to determine the acceptance. 
While less statistically powerful, this technique yields a value of Re(e'/e) which is consistent 
with the standard analysis based on the uncorrelated parts of the statistical and systematic 
errors. Finally, using 7c + tt~ data from 1996 (collected simultaneously with the 7r°7r° data) 
instead of from 1997 yields a value of Re(e'/e) which is consistent with the standard analysis, 
allowing a systematic error of 4 x 10~ 4 due to the 1996 Level 3 inefficiency. 

In conclusion, we have measured Re(e'/e) to be (28.0 ± 3.0 (stat) ± 2.8 (syst)) x 10~ 4 ; 
combining the errors in quadrature, Re(e'/e) = (28.0 ±4.1) x 10 -4 . This result definitively 
establishes the existence of CP violation in a decay process, agreeing better with the ear- 



lier measurement from NA31 than with E731 p2| |, and shows that a superweak interaction 



cannot be the sole source of CP violation in the K meson system. The average of the three 
measurements, (21.7 ± 3.0) x 10~ 4 , while at the high end of standard-model predictions, 
supports the notion of a nonzero phase in the CKM matrix. Further theoretical and exper- 
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imental advances are needed before one can say whether or not there are other sources of 
CP violation. 
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